Enhanced Gilbert Damping in Thin Ferromagnetic Films 
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The precession of the magnetization of a ferromagnet is shown to transfer spins into adjacent 
normal metal layers. This "pumping" of spins slows down the precession corresponding to an en- 
hanced Gilbert damping constant in the Landau-Lifshitz equation. The damping is expressed in 
terms of the scattering matrix of the ferromagnetic layer, which is accessible to model and first- 
principles calculations. Our estimates for permalloy thin films explain the trends observed in recent 
experiments. 
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The magnetization dynamics of a bulk ferromagnet is 
well described by the phenomenological Landau-Lifshitz- 
Gilbert (LLG) equation 0] 



dm 



-7m x H ff + am x 



dm 



(1) 



where m is the magnetization direction, 7 is the gy- 
romagnetic ratio, and H c ff is the effective magnetic 
field including the external, demagnetization, and crys- 
tal anisotropy fields. The second term on the right-hand 
side of Eq. (|l|) was first introduced by Gilbert jjj and the 
dimensionlcss coefficient a is called the Gilbert damp- 
ing constant. For a constant H c s and a — 0, m pro- 
cesses around the field vector with frequency u> = 7-ff e ff- 
When damping is switched on a > 0, the precession spi- 
rals down to a time independent magnetization along the 
field direction on a time scale of l/au>. Study of a in 
bulk metallic ferromagnets has drawn a significant in- 
terest over several decades. Notwithstanding the large 
body of both experimental and theoretical [|] work, 
the damping mechanism in bulk ferromagnets is not yet 
fully understood. 

The magnetization dynamics in thin magnetic films 
and microstructures is technologically relevant for, e.g., 
magnetic recording applications at high bit densities. Re- 
cent interest of the basic physics community in this topic 
is motivated by the spin-current induced magnetization 
switching in layered structures pi |L 0|. The Gilbert 
damping constant was found to be 0.04 < a < 0.22 for 
Cu-Co and Pt-Co || which is considerably larger than 
the bulk value «o ~ 0.005 in Co ^,^|. Previous attempts 
to explain the additional damping in magnetic multilayer 
systems involved an enhanced electron-magnon scatter- 
ing near the interface || and other mechanisms jl0| , both 
in equilibrium and in the presence of a spin-polarized cur- 
rent. 

In this Letter we propose a novel mechanism for the 
Gilbert damping in normal-metal-ferromagnet (N-F) hy- 
brids. According to Eq. ([I]), the precession of the magne- 



tization direction m is caused by the torque oc m x H e ff. 
This is physically equivalent to a volume injection of what 
we call a "spin current" . The damping occurs when the 
spin current is allowed to leak into a normal metal in con- 
tact with the ferromagnet. Our mechanism is thus the 
inverse of the spin-current induced magnetization switch- 
ing: A spin current can exert a finite torque on the fer- 
romagnetic order parameter, and, vice versa, a moving 
magnetization vector loses torque by emitting a spin cur- 
rent. In other words, the magnetization precession acts 
as a spin pump which transfers angular momentum from 
the ferromagnet into the normal metal. This effect can be 
mathematically formulated in terms of the dependence of 
the scattering matrix of a ferromagnetic layer attached 
to normal metal leads on the precession of m, analogous 
to the parametric charge pumping in nonmagnetic sys- 
tems JlT| . The damping contribution is found to obey 
the LLG phenomenology. Enhancement of the damping 
constant a' — a — ao can be expressed in terms of the 
scattering matrix at the Fermi energy of a ferromagnetic 
film in contact with normal metal reservoirs, which can 
be readily obtained by model or first-principles calcula- 
tions. Our numerical estimates of a' compare well with 
recent experimental results fl2|| . Earlier experiments re- 
ported in Ref. Jl3| can also be understood by our model 

B 

We consider a ferromagnetic film sandwiched between 
two paramagnetic layers as shown in Fig. |l|. Spin pump- 
ing is governed by the ferromagnetic film and the vicinity 
of the N-F interfaces. The normal metal layers are, there- 
fore, interpreted as reservoirs attached to non-magnetic 
leads. The quantity of interest is the 2x2 current matrix 
in spin space / = l/ c /2 — <r • I s e/ft, for the charge (I c ) 
and spin flow (I s ) from the magnetic film into adjacent 
normal metal leads, where I is the unit matrix and <r the 
vector of Pauli spin matrices. 

When no voltages are applied and the external field 
is constant, the charge current vanishes. Two contribu- 
tions to the spin current I s on cither side of the ferro- 
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FIG. 1: Ferromagnetic film (F) adjacent to two normal metal 
layers (N). The latter are viewed as reservoirs in common 
thermal equilibrium. The reflection and transmission ampli- 
tudes r and t' shown here govern the spin current pumped 
into the right lead. 



magnet may be distinguished, viz. IP um P an( j l^ ). jpump 
is the spin current pumped into the normal metal to be 
discussed below, whereas is the current which flows 
back into the fcrromagnet. The latter is driven by the 
accumulated spins in the normal metal and gives, e.g., 
rise to the spin-current induced magnetization switching 
j|, ||, q|. Here, we model the normal metal as an ideal sink 
for the spin current, such that a spin accumulation does 
not build up. This approximation is valid when the spins 
injected by IP um P decay and/or leave the interface suffi- 
ciently fast, e.g., when the dimensionless conductance of 
the N-F interface is smaller than h/T s {5 |Q. Here, r s f 
is the spin-flip relaxation time and S is the energy level 
spacing at the Fermi surface of a normal metal film with 
a thickness which is the smaller one of the geometrical 
film thickness and the spin-flip diffusion length. 

The current I{t) pumped by the precession of the mag- 
netization into the right and left paramagnetic reservoirs, 
connected to the ferromagnet by normal metal leads (R) 
and (L), may be calculated in an adiabatic approxima- 
tion since the period of precession 2tt /lu is typically much 
larger than the relaxation times of the electronic degrees 
of freedom of the system. The adiabatic charge-current 
response in nonmagnetic systems by a scattering matrix 
which evolves under a time-dependent system parameter 
X(t) has been derived in jn], Adopting Brouwer's 
notation | II , the generalization to the 2x2 matrix cur- 
rent (directed into the normal metal lead I — R or L) 
reads 



I(t) 



pump 



dn(l) dX{t) 



dX dt ' 
where the matrix emissivity into the lead I is 



dn{l) 



1 



\ " 



ds 



mn,W ~t 



dX Am ^ 

mnl 



H.c. 



(2) 



(3) 



and s is the 2x2 scattering matrix of the ferromag- 
netic insertion, m and n label the transverse modes 
at the Fermi energy in the normal metal leads and 
V = R,L. Spin-flip scattering in the contact is dis- 
regarded, s depends on the magnetization direction 
m of the fcrromagnet through the projection matrices 



= (i + & ■ m) /2 and = (i - <x ■ m) /2 (Ref. (l6|): 



Smn,ll> = S'mn.ll'^ + S mn,W Tn<3 S P in curren t pumped 

by the magnetization precession is obtained by identi- 
fying X(t) — <p(t), where (p is the azimuthal angle of 
the magnetization direction in the plane perpendicular 
to the precession axis. The resulting current is traceless, 
jpump _ _^ e jjj^fr . jpump^ j g ^ charge current indeed 

vanishes, and 



jpump 



h 

47T 



A r m x 



dm 



At 



dm 
~dt 



(4) 



where the interface parameters are 

A - i-Y^lM -r l | 2 + I/' T -t' 1 \ 2 \ 

— ey / j \\> mn mn\ T I'mn v mn\ J) \° I 

V r T ( r l )* i f'T U'i \*\ (6) 

, L mn\ ran) 1 mn\ mnJ J * \ w y 



^4, 
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Im 



inn 



mn 



' mnl 



is the reflection coefficient for spin-up 
[spin-down] electrons in the Ith lead and t'^ nn [i^„] is the 
transmission coefficient for spin-up [spin-down] electrons 
into the Zth lead. (See Fig. |l| for I = R.) Using unitarity 
of the scattering matrix for each spin direction, we can 
summarize Eqs. (||) and (||) by A r +iAi = g^~t^, where 

9^' = Emnt'Smn-CnKm)*} is the ( DC ) Conductance 

matrix TO @ and *ti = £nm The spin 

current H) trivially vanishes for the steady state, i.e., 
when dm/dt — 0, and for unpolarized contacts s] nn a , = 
s i , 

mn.ll ' 

Per revolution, the precession pumps an angular mo- 
mentum into an adjacent normal metal layer which is 
proportional to A r , in the direction of the (averaged) ap- 
plied magnetic field, and decaying in time. At first sight, 
it is astonishing that a pump can be operated by a sin- 
gle parameter varying in time, whereas the "peristaltic" 
pumping of a charge current requires at least two pe- 
riodic parameters |Tl] ]. However, there are actually two 
periodic parameters (out of phase by 7r/2) hidden behind 
<p{t), viz. the projections of the unit vector defined by ip 
in the plane perpendicular to the axis of precession. 

By conservation of angular momentum, the spin torque 
on the ferromagnet resulting from the spin pumping into 
the nonmagnetic leads gives an additional term to the 
LLG equation (|l|). After including this term, Eq. (Q) re- 
mains valid, but the gyromagnetic ratio and the damping 
constant are renormalizcd: 
I 
7 



-{l + g^+A^/AnM}. 

7o 



j_ 

To 



{ao+g L [4^+A^]/AirM}. 



(7) 
(8) 
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Here, g^ is the Lande factor and M is the total magnetic 
moment (in units of /is) of the ferromagnetic film; sub- 
script denotes the bulk values of 7 and a; superscripts 
(L) and (R) denote parameters evaluated on the left and 
right side of the F layer, respectively. Eqs. (^) and (||) 
are the central result of this paper. A r and Ai affect, 
e.g., ferromagnetic resonance experiments as a shift of 



the resonance magnetic field via A, 



A[ R \ whereas 



,(£) 



Ar increases the relative resonance linewidth. 



From now on we focus on ferromagnetic films which 
are thicker than the coherence length Af c = tt / (fc-j- — fcjj, 
where ku are the spin-up and spin-down Fermi wavevec- 
tors, i.e., thicker than a few monolayers in the case of 
transition metals. In this regime, spin-up and spin-down 
electrons transmitted or scattered from one N-F interface 
interfere incoherently at the other interface, tT J. vanishes 
and the mixing conductance g'^ is governed by the re- 
flection coefficients of the isolated N-F interfaces. 

Ai = Img'* vanishes for ballistic and diffusive con- 
tacts as well as nonmagnetic tunnel barriers ]lq| . First- 
principles calculations find very small Ai for Cu-Co and 
Fe-Cr It is, therefore, likely that Ai may be disre- 
garded in many systems. If Ai does vanish on both sides 
of the ferromagnetic film, it follows from Eqs. (0) and 
(|J) that the resonance frequency is not modified 7 = 70 
and the enhancement of the Gilbert damping is given by 
a' =g L [Ar L] + A [ r R) ]/AnM. 

The coefficient A r can be estimated by simple model 
calculations JHJ. For ballistic (point) contacts, A® — 
(1 +p)g with the polarization p = f«1t _ nil's 
and the average conductance g - 



(ff TT_ 5 U)/( ff TT +5 U) 
( 5 TT + g ll)/2. For dif- 
fusive N-F hybrids, A® = g^, the conductance of the 
normal metal part. A nonmagnetic tunneling barrier be- 
tween F and N suppresses the spin current exponentially. 
The magnetization precession of a magnetic insulator can 
also emit a spin current into a normal metal, since g^ 
does not necessarily vanish because the phase shifts of re- 
flected spin-up and spin-down electrons at the interface 
may differ Q . 

Let us now estimate the damping coefficient a' for 
thin films of permalloy (Ni 80 Fe2o, Py)j a magnetically 
very soft material of great technological importance. 
Mizukami et ai. ]T^ | measured the ferromagnetic reso- 
nance linewidth of N-Py-N sandwiches and discovered 
systematic trends in the damping parameter as a func- 
tion of Py layer thickness d for different normal met- 
als. The spin polarization of electrons emitted by Py has 
been measured to be p « 0.4 in point contacts ]l9| , the 
magnetization per atom is / w 1.2, and Lande factor- 
gL ~ 2.1 p5. The interface conductance of metallic 
interfaces with Fe or Co is of the order of 10 15 f2 _1 m~ 2 , 
with significant but not drastic dependences on interface 
morphology or material combination |20[ . This corre- 
sponds to roughly one conducting channel per interface 
atom. Assuming the Fermi surface of the normal metal is 
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FIG. 2: The lines show our theoretical result (9) with k =1.0, 
0.6, and 0.1; the data points are derived from the measure- 
ments [12] shown in two insets. Insets: Measured Gilbert 
damping constant a (lower inset) and the relative atomic 
magnetization f/fo (upper inset) in permalloy film of varied 
thickness dp y in a trilayer structure N-Py-N. 



isotropic, we arrive at the estimate a' rj l.l/d(A). The 
factor 1/d does not reflect an intrinsic effect; a reduced 
total magnetization is simply more sensitive to a given 
spin-current loss at the interface. Comparing with the in- 
trinsic a w 0.006 of permalloy (jl2], |2l| the spin-current 
induced damping becomes important for ferromagnetic 
layers with thickness d < 100 A. We can refine the es- 
timate by including the significant film-thickness depen- 
dence of the magnetization measured by the same group 
fll2fl. We, therefore, improve our above estimate as 



a'(d) 



1.1 fo 

K X -—-5— X 



d(A) f(d) 



(9) 



where fo and f(d) are the atomic magnetization of the 
permalloy bulk and films, k is an adjustable parameter 
representing the number of scattering channels in units 
of one channel per interface atom, which should be of the 
order of unity. 

The experimental results for the damping factor a and 
the relative magnetization f/fo for N-Py-N sandwiches 
with iV=Pt, Pd, Ta, and Cu are shown in the insets 
of Fig. ||. Our estimate (^) appears to well explain the 
dependence of a on the permalloy film thickness d (see 
Fig. |^) for reasonable values of n. First-principles calcu- 
lations are called for to test these values. 

The lack of a significant thickness dependence of damp- 
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ing parameter of the Cu-Py system requires additional 
attention. An opaque interface might be an explanation, 
but it appears more likely that due to long spin-flip re- 
laxation times in Cu, the 5 nm thick buffer layers in 
do not provide the ideal sink for the injected spins as 
assumed above. This means that a nonequilibrium spin 
accumulation on Cu opposes the pumped spin current 
and nullifies the additional damping when h/r s f5 is com- 
parable or smaller than the conductance g. For 5 nm Cu 
buffers, gS/h ~ 10 13 s" 1 , whereas l/r sf ~ 10 12 s" 1 @. 
It follows that Cu is indeed a poor sink for the injected 
spins and the Gilbert damping constant is not enhanced. 
On the other hand, Pt, Ta, and Pd are considerably heav- 
ier than Cu and, since l/r S f scales as Z 4 |2^| , where Z 
is the atomic number, have much larger spin-relaxation 
rates and our arguments hold. 

A physical picture of the effect of magnetization preces- 
sion in layered systems has been proposed earlier by Hur- 
dequint et al. |24f| in order to explain ferromagnetic and 
conduction electron spin resonance experiments. These 
authors realized that the precessing magnetization is a 
source of a nonequilibrium spin accumulation which dif- 
fuses out of the N-F interfaces into the adjacent normal 
metal layers where it can dissipate by spin-flip processes. 
Enhanced Gilbert damping in thin ferromagnetic films 
in contact with normal metal has also been discussed by 
Berger for a ballistic N-F interface in a spin-valve 
configuration. His expression for the damping coefficient 
[Eq. (20) in Ref. Bj] scales like ours as a function of layer 
thickness, but differs as a function of material parame- 
ters. E.g., in contrast to our result, Berger's expression 
does not vanish with vanishing exchange splitting. 

In conclusion, we demonstrated that the Gilbert damp- 
ing constant is enhanced in thin magnetic films with nor- 
mal metal buffer layers by a spin-pump effect through the 
N-F contact. The damping is significant for transition 
metal films thinner than about 10 nm. Recent experi- 
ments on permalloy films |fl2f are well explained. 

We are grateful to L. Berger, B. I. Halpcrin, B. Hcin- 
rich, D. Huertas-Hernando, and Yu. V. Nazarov for stim- 
ulating discussions. This work was supported in part by 
the DARPA Award No. MDA 972-01-1-0024, the Nor- 
wegian Research Council, the NEDO International Joint 
Research Grant Program "Nano-magnetoelectronics" , 
NSF Grant No. DMR 99-81283, and the Schlumberger 
Foundation. 

Note added. — After submission of this paper, we be- 
came aware of new exciting experimental results |2j| on 
enhanced Gilbert damping in ultrathin iron films which 
support the spin-pumping mechanism proposed in this 
Letter |B3. 
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